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We present an experimental study of a droplet interacting with an ultrasonic surface acoustic 
wave (SAW). Depending on the amplitude of the wave, the drop can either experience an internal 
flow with its contact-line pinned, or (at higher amplitude) move along the direction of the wave 
also with internal flow. Both situations appear together with oscillations of the drop free-surface. 
The physical origins of the internal mixing flow as well as the drop displacement and surface waves 
are still not well understood. In order to give insights of the underlying physics involved in these 
phenomena, we carried out an experimental and numerical study. The results suggest that the 
surface deformation of the drop can be related as a combination between acoustic streaming effect 
and radiation pressure inside the drop. 

PACS numbers: 47.55.D- ; *43.25.Nm ; 68.08.Bc ; *43.25.Qp 



I. INTRODUCTION 
A. General Issues 

For discrete microfluidics in lab-on-chip devices, it is 
often necessary to handle small amount of liquids in 
the form of droplets. These elementary operations are 
the linear displacement, the splitting or the merging of 
droplets, or the production of a constant flow inside a 
drop. At the scale of a micro- litre drop or smaller, the 
main parameters governing the shape and flow dynam- 
ics, are the dynamic viscosity \i and the surface tension 
a. Non-linear effects originating from inertia are mostly 
negligible, so that the momentum conservation equation 
for the fluid reduces to the time-dependent Stokes equa- 
tion. Consequently, it is hard to induce (chaotic) mix- 
ing and continuous re-suspension of particles in a drop 
as the equation of flow motion is time-reversible. Fur- 
thermore, during these operations, the drop is often in 
contact with ambient atmosphere and evaporation can 
significantly occur at the typical time-scales of the op- 
erations. It is known that when the liquid contains col- 
loidal particles or macromolecules, evaporation leads to 
a outwards flow in the drop, which in turn leads to a 
particle clustering in the vicinity of the contact-lineplQ. 
The 'coffee rings' originate from this phenomenon, which 
in practical situation can strongly limit bio-particles de- 
tection in lab-on-chip devices (see, e.g., Q) by forming 
deposits on DNA microarrays [4(. 

Another difficulty is to displace such small droplets. 
As their typical size is smaller than the capillary length 

= v 1 o 'i 1 pg, of the order of 1 to 2 mm for most liquids 
(p being the liquid density), the volume forces are gen- 
erally overcome by the retention force, Fr, acting at the 
contact-line of the substrate. This force scales with the 
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FIG. 1: (a) Successive snapshots of a water drop displaced 
by SAW, showing a periodic creeping and jumping flow. The 
arrow shows the direction of the displacement, (b) Particle 
tracking inside a steady drop subjected to SAW waves, ob- 
tained by time- averaging, (c) Asymmetric shape of a drop 
moving under the action of SAW (W/G mixture 2), without 
oscillations. For the three images, the SAW propagates from 
left to right. 



contact-line perimeter 



Fr = 7rR<j(cos(6 r ) — cos(O a )) 



(i) 
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Its strength is related to the contact angle hysteresis 
a — r , where a and r are respectively the advanc- 
ing and receding contact angles. Real surfaces combin- 
ing micro- or nano-scale roughness and chemical hetero- 
geneities, can have significant hysteresis [6|. Further- 
more, the aforementioned accumulation of particles near 
the contact-line increases this hysteresis even more. The 
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coffee-ring effect can be hindered by Marangoni forces if 
one uses a very volatile solvent [7], but with casual liq- 
uids and substrates used in biology, volume forces are too 
weak to re-suspend continuously the particles or to move 
the drop in a controlled fashion. 

One of the challenging issue, is to provoke the contin- 
uous unpinning of the contact-line of the drop to both 
prevent particle accumulation and to obliterate hystere- 
sis. Various authors evidenced a similar unpinning effect 
with low- frequency mechanical vibrations [8l4l3j , but it is 
required that the amplitude of vibrations be of the order 
of a millimetre to obtain the desired effect. Such vibra- 
tions can be unsupported by many fragile surrounding 
devices. Furthermore, as the pinning forces are often due 
to micron-sized roughness, it is intuitively more efficient 
to induce an unpinning flow from waves of shorter wave- 
length such as those produced by ultrasonic transducers, 
A ranging from 1 to 300 microns. 



B. Qualitative description of the phenomena 

Recently, it has been found that the use of ultrasonic 
surface acoustic - so-called 'Rayleigh' - waves (SAW) 
transducers on piezoelectric substrates enables to both 
agitate inner fluid and to displace a droplet along the di- 
rection of the wave propagation (TUTSI. SAW frequency 
generally ranges from 5 to 150 MHz, and depending on 
the operating frequency the dynamics of the droplet can 
show a host of different behaviours. At low frequencies 
- about 20 MHz, a 'creeping and jumping' motion of the 
drop generally occurs, especially for low-viscosity liquids 
[17j, see Fig. [l]-a. Due to a non-linear coupling between 
the Rayleigh acoustic wave and the flow, the liquid inside 
is highly stirred in a rotating flow motion (see Fig. [l]-b), 
while the interface shows capillary waves. Although sev- 
eral experimental studies demonstrated the reliability of 
this technique, little is known about the detailed mecha- 
nisms that create both the internal flow and the droplet 
deformations. It is generally admitted that the SAW is 
radiated inside the droplet at a refraction angle <p r , anal- 
ogous to Snellius-Descartes law in optics: 

sin(0 r ) = ^ (2) 

where c\ and c s denote respectively the speed of sound in 
the liquid and the solid, and in turn generates acous- 
tic streaming. The angle (j) r is around 25° for water 
and usual piezo-electric substrates. Regarding droplet 
dynamics by SAW, a fluid dynamics analysis has been 
recently attempted [15j , but results are mostly related to 
atomisation occurring at high acoustic power. At high 
frequencies - around 100 MHz, a noticeable contribution 
on the understanding of drop deformation and displace- 
ment, is available in [H| (however at this frequency, no 
oscillations have been observed). It shows that the flow 
originates from a carrying force focused in a narrow re- 



gion inside the drop, starting at the drop edge hit by the 
wave and directed towards the refraction angle. 

In the general case, this force is partly due to non-linear 
acoustic streaming [19|, [20|. The fluid trajectories are 
described qualitatively by the particle tracking snapshot 
like that of Fig. QJb. In the present frequency range - 
around 20 MHz, the liquid- air interface is also affected 
by SAW in three different ways: 

- First, the droplet 's left-right symmetry is broken 
(see Fig. [IJc) and the asymmetry is more pronounced 
for larger acoustic power. The asymmetry is responsi- 
ble for the drop displacement as it makes the front and 
rear contact-angles becoming respectively larger than a 
and smaller than r [6[. However, the relative contribu- 
tions of acoustic streaming and radiation pressure on this 
asymmetry with respect to the frequency regime, remain 
unclear. 

- Second, the periodic creeping-jumping motion (see 
Fig. [T]-a) has a well-defined frequency of the order of 50 
to 200 Hz for usual droplet size, and comes with a contin- 
uous unpinning of the contact-line. During the deforma- 
tion cycle, the droplet is strongly stretched and flattened. 

- Finally, besides the global deformation, the droplet 
surface exhibits small deformations ('trembling') at ap- 
proximately regular frequency of a few thousands Hertz. 
However, this is hardly visible on the pictures presented 
here. 

In any case, the periodic deformations of the inter- 
face cannot be attributed to a steady flow. Instead, the 
acoustic radiation pressure which affects the shape of the 
liquid- vapour interface, can be invoked for a possible ori- 
gin of such deformations. Therefore, the periodic motion 
could be associated with dynamical eigenmodes of the 
droplet, where capillarity tends to restore a minimal sur- 
face energy. 

The deep understanding of the involved phenomena 
is necessary in order to set-up a droplet handling lab- 
on-chip with minimal acoustic power towards an optimi- 
sation of the system. Especially interesting is the use 
of several combined transducers and short-time pulses, 
in order to reduce heat (continuous appliance of SAW 
would cause much heat production) and to protect fragile 
bio-particles. The current available literature being in- 
sufficient for quantitative predictions, the present study 
aims to carry out quantitative measurements of the dis- 
placement and deformation of the drop. We present here 
such a quantitative study by varying parameters like the 
droplet volume, the liquid viscosity and the amplitude of 
the acoustic wave. These results constitute a first step to- 
wards the understanding of the detailed coupled acoustic 
and hydrodynamic mechanisms. 



II. EXPERIMENTAL SETUP 

The set-up is depicted in Fig. [2j We used a sub- 
strate with piezoelectric properties (Lithium Nyobate, 
LiNbOs), in order to generate powerful SAW. The trans- 
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Lithium Niobate Substrate Fingers Drop 
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Liquid Kin. viscos. v Surf, tension 7 Density p 

(mm 2 /s) (N/m) (g/cm 3 ) 

Water 1.00 0.072 1.00 

W/G. Mix 1 5.60 0.067 1.126 

W/G. Mix 2 11.50 0.066 1.161 

TABLE I: Physical properties of liquids 



varied from 50 to 70, increasing the viscosity by a factor 
of up to 11, whereas the surface tension stayed close to 
that of water. 



FIG. 2: Scheme of the SAW interdigitated transducer that 
actuates with a droplet deposited on the substrate. The wave 
propagates from left to right. Below: The details of the net- 
work of fingers that compose the IDT. 



ducer is an inter-digitated transducer (IDT) which gen- 
erates transverse acoustic waves propagating along the 
surface. The interconnected fingers of the IDT are made 
of Titane covered by Gold, and designed by a lithographic 
technique detailed elsewhere [14[. We applied the peri- 
odic sinusoidal voltage with a high-frequency generator 
(IFR 2023 A), amplified with a home-made amplifier. A 
SAW is generated providing that the voltage frequency is 
compatible with the space between each track of the IDT. 
Considering that the space between fingers a/2 is 43.75 
/im, giving a wavelength A = 2a = 175 /im, and that the 
sound velocity in LiNbC>3 is 3485 m.s -1 for transverse 
waves, the value of the frequency fo has to be around 
19.5 MHz. In practice, we found the best actuation to 
droplets of any size for fo = 20.375 MHz and we kept this 
value for all experiments. The properties of the SAW en- 
sure that the amplitude is not attenuated in its direction 
of propagation, and that it is localised near the surface. 

In order to reduce the friction forces, and to study 
the simpler situation of partial wetting, we treated the 
surface with hydrophobic coating (monolayer of OTS). 
The contact angles are: # a =105° and # r =95°. This weak 
hysteresis helps to conduct experiments without the ad- 
ditional complexity of pearling at the trailing edge of 
a moving drop [22| that would have occurred otherwise. 
The droplet dynamics is acquired with a high-speed cam- 
era (Photron SA3) with a zoom lens completed by exten- 
sion rings. This allows a maximal magnification of about 
8 /im/pixel. The acquisition rate ranged from 2000 to 
5000 im/s. To avoid evaporation, special care is devoted 
to use cold source of light (Schott KL2500) and to acquire 
pictures no later than a few seconds after the droplet has 
been put on the substrate. 

The liquids used are deionised water and wa- 
ter /glycerol mixtures (W/G), which physical properties 
at 20° C are given in Table 1. Glycerol percent weight was 



The measurement of the normal displacement at the 
surface of the substrate gives the amplitude of the SAW. 
We used a technique of laser heterodyne interferometer 
of type Mach-Zender (SH130, B.M. Industries). In brief, 
the interferometer measures the Doppler shift between 
a reference beam and a secondary beam modulated at 
Im = 70 MHz by an acousto-optic Bragg cell. The sec- 
ondary beam is reflected onto the surface and hence, is 
modulated a second time by Doppler effect. The acoustic 
displacement d is then deduced from these measurements 
by extracting the relative amplitude between the peaks 
at /m and /m ± fo i n the power spectrum. The displace- 
ment d ranges from 0.3 to 2 nm for the usual experimental 
conditions. 



III. RESULTS 

For each experiment, the dynamics of the drop was ex- 
tracted using the public domain image processing soft- 
ware Image J (http://rsbweb.nih.gov/ij/). From a se- 
quence of about 2000 images, we built spatio-temporal 
diagrams that allowed us to measure both the speed and 
the frequency of oscillations: we extracted the grey lev- 
els along a horizontal line cutting the base of the droplet, 
and we pilled the lines on top of each others (time axis 
is horizontal after a rotation of tt/2). A typical motion 
of drop is pictured in Fig. [3] after binary treatment. The 
total time is about 250 ms. This diagram is then skele- 
tonised to get the position versus time. The velocity U is 
the slope and the frequency / is extracted from a Fourier 
transform. It is noticeable that the drop does not start 
to move or oscillate immediately: it waits about 200 ms 
before it starts. During this relatively long transient, the 
oscillations progressively increase, as does the velocity. 



A. Droplet displacement 

We first present the velocity measurements, obtained 
by varying the acoustic amplitude d at constant volume 
(Fig. [4]-a) and by varying the droplet volume V at con- 
stant displacement d (Fig. HJ-b), for water. As stated 
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FIG. 3: Typical spatio-temporal diagram for a moving, creep- 
ing and jumping droplet like in Fig. [I}a. 



earlier, the drop displacement results from a balance be- 
tween driving forces - due to both the streaming and 
radiation pressure which induce a left / right asymmetry - 
and dissipation occurring near the contact-line. 

Influence of acoustic displacement - As expected, 
the larger d is, the faster the droplet moves. Also it 
turns out that d has to be larger than a certain threshold 
around 0.5 nm in order for the droplet to move. Below 
this threshold, the droplet stays at the same location 
but an internal flow motion is observed: this situation 
is illustrated in Fig. [T]-b. The influence of d is strongly 
non-linear: after a sharp increase just above threshold, 
U saturates at a volume-dependent value. 

Influence of volume - Whatever the amplitude of the 
acoustic wave d, the influence of volume shows a maximal 
velocity at V about 5 /A (see Fig.|H-b). A tentative expla- 
nation of this behaviour is developed in the Discussion 
part. 

Influence of viscosity - To have a more complete 
understanding of the coupled acoustic and hydrodynamic 
effects, we used liquids of various viscosity (see Table 1). 
Figure [5] shows velocity versus volume, for four liquids 
(Water and W/G Mixtures 1 and 2). The increase of 
viscosity leads to a decrease of the velocity as expected 
for moving droplets. Indeed it is known that viscosity is 
involved in the dissipation near the contact-line [6] and, 
hence, an equal driving force displaces viscous drops at 
lower speed. It is natural to consider the dimensionless 
velocity, the capillary number 



Ca=^ 
a 



(3) 



which is plotted in the insert of Fig. [5] 

It turns out that the data for capillary numbers do not 
collapse well together, although the order of magnitude 
is around lxlO -3 for most measurements. This can be 
interpreted in the following way: the driving force is due 
(at least partially) to the acoustic streaming which in- 
volves bulk dissipation [20]. Hence, streaming effects are 
promoted by liquid viscosity. The increase of viscosity 



(a) 



(b) 



o V = 2 |il 
• V = 5ul 

□ V= 15 mJ 



□ □ 



• ft 

□ 



an 



o 

• 

o 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

d (nm) 





100 




90 




80 




70 




60 


CO 




E 






50 


z> 


40 




30 




20 




10 








□ 



\- CD 
O 



Oo° o 

o o 



o d = 0.56 nm 
• d = 1 .17 nm 

□ d = 1.86 nm 



10 



15 20 

voii) 



25 



30 



FIG. 4: (a) Drop velocity versus acoustic displacement d, for 
3 different volumes (water), (b) Drop velocity versus volume, 
for 3 values of d. 



increases the contact-line dissipation [6], but it should 
also increase the driving force. Therefore, it is consistent 
that the values of Ca be slightly larger for more viscous 
liquids. 



B. Creeping-jumping oscillatory motion 

Influence of viscosity - First, it is noticeable that 
oscillations are more and more damped as the viscos- 
ity is increased by adding glycerol in the mixture. For 
instance with mixture 1, oscillations have very small am- 
plitude and are barely measurable, whereas for mixture 
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FIG. 5: Droplet velocity versus volume, for three different 
values of viscosity, d = 0.99 nm. Inserted: the same measure- 
ments with the dimensionless velocity Ca. 
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FIG. 6: Frequency of the creeping-jumping oscillations versus 
volume, for various acoustic displacement d. The dashed blue 
curve is the best fit / ~ V _1//2 , and the black plain curve is 
the prediction by Rayleigh 's theory. 



fn 



n(n - l)(n + 2)7 
3npV 



1/2 



(4) 



in which f n denotes the resonance frequency of the n th 
mode of oscillation. 

The case n=l corresponds to pure translation and n=2 
corresponds to the lowest mode in which the drop is elon- 
gated with an elliptical shape. The morphology of the 
oscillations suggests a mode n=2 (see Fig. [I]-a). To the 
best of our knowledge, the case of sessile drops has not 
been treated fully analytically. The pinning force acting 
at the contact-line strongly constrains the structure of 
the capillary wave, making its computation more com- 
plex. However, Noblin et al. [TTJ] showed that for mod- 
erate hysteresis and under strong enough vibrations, the 
contact-line is constantly unpinned and the constraint 
due to retention force at the contact-line (eq. jTJ) is re- 
leased 

110, EMU- 

This is indeed what is observed in 
sequences like Fig. [l]-a. Hence, we compared out results 
to what is predicted by eq. (jlj, also plotted in Fig. [6l 
It turns out that Rayleigh-Lamb's theory slightly over- 
estimates the results although the order of magnitude is 
correctly predicted. This can be interpreted by the fact 
that SAW tend to induce large amplitude droplet oscil- 
lations which are out of the frame of Rayleigh's theory. 
Intuitively, large amplitude oscillations are of larger pe- 
riod than weak amplitude ones, and this is confirmed by 
the results of Fig. El larger d generates free-surface oscil- 
lations of larger amplitude and of smaller frequencies. 

Amplitude - In Fig. [7J the relative amplitude of os- 
cillations - i.e. the amplitude divided by the mean base 
radius of the drop - is plotted versus volume, for differ- 
ent acoustic displacements d. This amplitude is extracted 
from spatiotemporal diagrams like that in Fig. [3] and cor- 
responds to the maximal lateral deformation of the base 
of the drop. It turns out that the amplitude is maximal 
for a d-dependent value of volume. These results suggests 
that oscillations result from a balance between the wave- 
induced pressure reaching the interface and capillarity, 
as developed in the Discussion. 



2 oscillations completely vanished. Hence, the results of 
oscillations reported thereafter are for water as this is the 
only situation where they clearly appear like in Fig. [IJa. 

Frequency - In Fig. [6j / is plotted versus volume for 
different acoustic displacement d. 

The results clearly show a sharp decrease of / with V, 
in the range V = 0.5 to 5 fA. A power-law fit leads to 
an exponent close to one half: / ~ V~ x l 2 . This is rem- 
iniscent to what is predicted by the theory of Rayleigh- 
Lamb [2l|, for inertial-capillary modes of vibrations in a 
drop. Rayleigh's seminal calculation is valid for a spher- 



IV. DISCUSSION - NUMERICAL RESULTS 

The measurements presented here show several evi- 
dences that the dynamics of the drop is due to an in- 
terplay between acoustic streaming and radiation pres- 
sure. For water droplets at the tested frequency, the SAW 
induces both a left-right symmetry break-up leading to 
motion and the creeping-jumping oscillations (the high- 
frequency 'trembling' oscillations, mentioned earlier, will 
be treated in a further study). It is still unclear what 
the respective weight of streaming and radiation pres- 
sure are on the drop left/right asymmetry. Schindler et 
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FIG. 7: Amplitude of the creeping-jumping oscillations of wa- 
ter droplets versus volume, for various acoustic displacements. 
The curves are guidelines for the eye. 



FIG. 8: The different regimes of the relative influence of ra- 
diation pressure versus frequency. The plain curve stands for 
the attenuation length L in pure water (cz =1476 m/s and 
/x b =2.8xlCT 3 Pa.s) given by eq. (©. 



al. [18[ showed that at high frequencies, acoustic stream- 
ing can lead to a conservative force potential in the drop 
towards the direction of propagation, headed to the top 
and inducing a deformation like in Fig.QJc. However, the 
streaming flow cannot explain the creeping-jumping os- 
cillations in which the free-surface of the drop is pushed 
upwards. 

The acoustic streaming inside the droplet comes from 
dissipation due to liquid shear (/x) and bulk vis- 
cosities. The internal streaming flow is more intense for 
larger viscosities and frequency [20] . On the contrary, the 
radiation pressure is directly proportional to the energy 
which reaches the free-surface [20|. Consequently, the 
larger the viscous dissipation is, the smaller this energy 
is. The relative importance of both effects is strongly de- 
pendent on the structure of the acoustic wave diffracted 
in the drop, and on how much the wave is attenuated. 

It is relevant to compare the typical drop size with the 
length of attenuation of the acoustic wave which, in a 
fluid, is given by (20| : 

L 2 ? c f ( 1 ^ ( 5 ) 

(27T/0) 2 If/i + wJ [ ' 

This length is plotted versus /o in Fig. [8j with respect 
to the drop size. The length L is equal or lower than 
the size of droplets for a frequency of about 125 MHz 
or above (see Fig. [8j) and defines a range of frequency 
where radiation pressure is negligible. This is far beyond 
the values of fo used here and, therefore, both streaming 
and radiation pressure are effective in our study. 

To confirm this trend, at least qualitatively, we carried 
out numerical simulations with a finite-elements software 
(COMSOL Multiphysics 3.4). The model uses a drop of 



half-spherical fixed shape (contact angle = 90 °, radius 
= 2 mm) subjected to a SAW. The General Form Modes 
of COMSOL Multiphysics with two domains has been 
used, one for the piezoelectric substrate and one for the 
water droplet. In the water droplet, the linear equation 
of acoustic in a thermo- viscous fluid has been considered. 
The boundary conditions used at the drop/ambient air 
interface are stress- free conditions. As we are interested 
in the linear radiation of the SAW propagating at the 
substrate surface into the water, the continuity of the 
normal stresses and displacements has been assumed be- 
tween the piezoelectric substrate and the liquid [23| . This 
last point is justified by the fact that in our experiences 
the amplitude of SAW is always small in comparison to 
the acoustic wavelength. 

Figure [9] shows the pressure distribution (a) in a water 
drop for /o=20 MHz and in liquids where the attenua- 
tion is (b) 25 times larger and (c) 100 times larger, for 
the same fo. When weakly attenuated, as in water for a 
frequency of 20 MHz (Fig.[9]-a), the wave structure is gen- 
erally chaotic due to multiple reflections as the drop acts 
like a chaotic cavity. If the radius is slightly modified, 
the wave structure is completely different except for the 
coherent part in the direct wave. To extract this coher- 
ent part we computed the wave structure for 15 droplets 
of radius ranging from 2 to 2.35 mm, and we averaged 
the contribution of each case. This calculation proce- 
dure, where stationary wave pattern are determined in 
non-deformable drops of increasing size, is justified by 
the fact that the time required for the wave to propagate 
across the 2 mm drop is about 2 /is, which is much lower 
than the observed characteristic deformation time of the 
drop. 
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FIG. 9: Numerical simulations of a SAW inside an non- 
deformable half-spherical drop, fo =20.375 MHz and d— 2 nm. 
Colour levels stand for pressure, (a) Weak attenuation in the 
case of water, (b) Intermediate attenuation, 25 times larger 
than (a), (c) Strong attenuation, 100 times larger than (a). 



The resulting data exhibits the coherent part of the 
wave, as shown in FiglTOl In this case, the wave reaches 
the free-surface whereas in the two other cases (Fig. [9]-b 
and c) it is strongly attenuated and hardly reaches the 
free-surface. It is clear that radiation pressure will only 
operate in the first case. The second case is streaming- 
dominated and should not show global free-surface oscil- 
lations. The action of radiation pressure is reminiscent to 




FIG. 10: Pressure due to the coherent direct wave in the situ- 
ation of Fig. [9}a, with chaotic multiple reflections. Obtained 
by averaging over several runs, see text for details. 



another study [24] , although for much smaller frequencies 
and a different set-up, where the free-surface of a liquid 
layer is deformed by acoustic vibrations of a plate. 

The coherent part of the wave is highly left-right asym- 
metric and the associated radiation pressure therefore 
contributes to the asymmetry of the droplet. On the 
contrary the chaotic part of the wave induces a radia- 
tion pressure homogeneously distributed on the surface 
of the wave. Therefore, only the direct coherent wave 
contributes to the displacement of the droplet whereas 
the chaotic part should be the driving effect of the free- 
surface oscillations. The computations showed that the 
chaotic part vanishes for smaller fo than the coherent 
direct part. Indeed, for the existence of multiple reflec- 
tions it is required that the length L be several times 
larger than the drop size. Therefore, we defined two other 
regimes: a 'multiple reflections' domain at lower /o, and 
an intermediate regime where only the coherent part of 
the direct wave contributes to the radiation pressure. 

The use of liquids more viscous than water have similar 
effects on the relative contribution of radiation pressure 
as the prescription of a higher frequency: it decreases L. 
In the cases displayed in Fig. [9]-b and Fig. [9]-c, the corre- 
sponding frequencies would be 100 MHz and 200 MHz, 
respectively. The use of W/G mixtures increases both \i 
and /ifc. The ratio between /i& and \i is 2.8 for water, and 
it should be quite close to this value for W/G mixtures 
of low viscosity like ours. It is sensible to interpret the 
absence of oscillations for W/G droplets by a decrease of 
the length L: as the wave is more strongly attenuated 
when arriving at the droplet free-surface, the radiation 
pressure is weaker than for water drops. In Fig. [8j the 
curve of L vs. fo for W/G mixtures would be just trans- 
lated downwards. 

To observe the creeping-jumping oscillations, the ra- 
diation pressure has to be strong enough to cause the 
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unpinning of the droplet contact-line (otherwise, as the 
drop volume is constant, the height of the drop cannot 
vary significantly). Our measurements and computations 
suggest that the chaotic part of the wave is required to 
obtain such pressure: this condition is fulfilled only for 
frequencies included in the lowest range domain and for 
low-viscous liquids. This is to be related to the results of 
Fig. [6] for the selection of the frequency: once the radia- 
tion pressure pushes the top of the drop and shrinks its 
base surface, the contribution of this pressure decreases 
as the height of the drop increases. Therefore, the drop 
retracts under the action of capillary forces, leading to 
a frequency comparable to that of the first eigenmode 
of the drop driven by an inertial-capillary balance. The 
maxima of amplitude found for a <i-dependent intermedi- 
ate volume (Fig. [7j) is consistent with this interpretation: 
surface tension hinders large oscillations for small drops, 
as capillary pressure, scaling with interface curvature, is 
bigger; whereas large drops are subjected to lesser radia- 
tion pressure because the acoustic wave is more damped 
before reaching the interface. 

Finally, these calculations suggest a possible explana- 
tion for the maximal drop velocity for water, obtained 
for intermediate drop volume (see Fig. |4]). If the drop 
is too small, only a part of the SAW is transferred to 
the droplet whereas a significant part of the acoustic en- 
ergy goes on propagating in the solid after the drop [25| . 
On the contrary, drops of larger size collect most of the 
acoustic energy for their displacement and deformation. 
Bigger drops collect about the same amount of energy as 
intermediate drops, but they will be slowed down simply 
because they have more inertia. The drop optimal size 
is dictated by the attenuation length in the solid, in the 
X direction. This attenuation length is no more linked 



to the thermo-viscous absorption in the liquid, but to 
the imaginary part ki of the Rayleigh wave wavenumber, 
which is of leaky type when propagating in a substrate in 
contact with a liquid. For a Lithium Niobate substrate 
loaded by water ki c± 2.3 10-5 x/o Np/m [23|, |26|, more 
than two orders of magnitude higher than the attenua- 
tion coefficient at fo=20 MHz. This corresponds to the 
fact that the Rayleigh wave efficiently looses its energy by 
radiating in the liquid. The extraction of the amplitude 
at the solid-liquid interface gives an attenuation length of 
about 2.13 mm, corresponding to a drop volume of about 
5 jil. This may explain the maxima of velocity observed 
in Fig. |H 



V. CONCLUSION 

We have presented here a quantitative study of droplet 
dynamics, when actuated by a SAW. Most of previous 
studies of SAW ignored the effect of radiation pressure 
as they were carried out at higher frequencies. In our 
range of frequency and liquid viscosity, a combination 
between streaming and radiation pressure effects was ob- 
served and their respective contributions on the observed 
dynamics was studied. The free-surface oscillations are 
mostly due to radiation pressure, whereas internal flow is 
due to streaming. This is evidenced by both experiments 
at different viscosity, and computations. Contributions 
of both phenomena are of interest, for an optimal micro- 
fluidics mixing efficiency. Future studies will focus on the 
influence of frequency in order to determine the relative 
weight of both mechanisms in the drop's asymmetry. 
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